Measuring the extragalactic background light from very high energy gamma-ray 

observations of blazars 
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The extragalactic background light (EBL) contains important information about stellar and 
galaxy evolution. It leaves imprint on the very high energy 7-ray spectra from sources at cos- 
mological distances due to the process of pair production. In this work we propose to measure the 
EBL directly by extracting the collective attenuation effects in a number of 7-ray sources at different 
redshifts. Using a Markov Chain Monte Carlo fitting method, the EBL intensities and the intrinsic 
spectral parameters of 7-ray sources are derived simultaneously. No prior shape of EBL is assumed 
in the fit. With this method, we can for the first time to derive the spectral shape of the EBL 
model-independently. Our result shows the expected features predicted by the present EBL models 
and thus support the understanding of the EBL origin. 

PACS numbers: 98.38.Mz,95.85.Ry 



The extragalactic background light (EBL) is the diffuse 
radiation from ultraviolet to far infrared wavelengths, 
spread isotropically in the universe (for a review of EBL, 
see [H, The EBL originates from the radiative en- 
ergy releases of all the stars, other extragalactic sources 
and diffuse emissions since the epoch of recombination. 
Therefore its intensity and spectral shape hold crucial 
information about the formation and evolution of stellar 
objects and galaxies throughout the cosmic history. The 
EBL is one of the fundamental quantities in cosmology. 

Direct measurement of EBL is, however, very difficult 
due to the contamination of the foreground emission from 
the solar system zodiacal light and the Galactic stellar 
and interstellar emissions . Technically, it also requires 
the absolute calibration of the instruments, and the un- 
derstanding all measurement uncertainties. Given the 
difficulties, direct measurements provide just lower and 
upper limits of EBL intensity. A strict lower limit on the 
EBL intensity is provided by the integrated light from 
resolved galaxies, e.g. in optical by the Hubble Space 
Telescope 3 and in infrared by the Spitzer telescope Q . 
The upper limit can be derived from the absolute mea- 
surement of EBL within its errors [l[ . The allowed range 
is shown by the shadow region in the following figures. 

Extreme efforts had been paid to calculate the EBL 
intensity The models generally include two dis- 

tinctive processes, that the UV and optical component 
of EBL is the integral of starlight over all epochs and 
the infrared component is due to the process of absorp- 



tion and re-emission of starlight by the interstellar dust. 
These models agree on the overall EBL shape, including 
two maxima at ^ Ifim by starlight and at ^ 100/^m by 
dust. However, since the detailed EBL model depends 
on many factors, such as the star formation history, the 
stellar initial mass function, the evolution of metallicity, 
the energy released by AGN, the size distribution and 
composition of dust grains, different models keep a large 
diversity. 

There is another indirect but effective way to study 
the EBL by observation of very high energy (VHE) 7- 
rays. The VHE 7-rays from extragalactic sources are 
attenuated by the process of electron/positron pair pro- 
duction, 7vHE + 7ERT, e+e~, when propagating to 



the Earth (e.g., ISMlSj). With the rapid development 



of ground based 7-ray imagining atmospheric Cerenkov 
telescopes (lACT), quite a few VHE 7-ray sources from 
cosmological distances have been detected, most of which 
are blazars, a subgroup of active galactic nuclei (AGN), 
with relativistic jet pointing towards the observer. With 
assumption of the intrinsic blazar spectra we can set an 
upper limit of the EBL intensity by cornparing the ob- 
served spectra with the intrinsic spectra [15[ . The obser- 
vations of blazars H 2356-309 and lES 1101-232 at red- 
shifts z = 0.165 and z = 0.186 respectively by HESS has 
set a strong upper limit of EBL, close to the lower limit 
set by galaxy counts, at the near infrared wavelength 
di]. The MAGIC observation of 3C 279 at z = 0.536 set 
upper limit at the optical band [l3]. In [l^ Mazin and 



2 



Raue gave a comprehensive study of EBL based on eleven 
blazars over a redshift range from 0.03 — 0.18. They ex- 
plored a large number of hypothetical EBL scenarios and 
set robust constraints on EBL over a wide wave-length 
range. With the Fermi observation of blazar spectra at 
GeV to ~ 100 GeV more stringent constraints on EBL 
are recently given by [iB-H. Those studies seem to in- 
dicate that the Universe is more transparent than we had 
expected. 

The power of this method to study EBL is limited due 
to the fact that the intrinsic spectrum of each blazar is 
unknown. Therefore it is hard to disentangle the absorp- 
tion effect by EBL from the intrinsic emission nature for a 
specific observation. The usual practice in the literature 
is to reconstruct the blazar intrinsic spectrum from the 
observation by first assuming an EBL model. The EBL 
model is rejected if it results in an unphysical intrinsic 
spectrum, for example, the reconstructed intrinsic spec- 
trum follows a power law with an extremely hard spectral 
slope or even shows an exponential rise at the high en- 
ergy end. Recently with large sample of 7-ray blazars, 
the EBL intensities were derived through a likelihood fit 
with given spectral template of the EBL 25 , 2^ . 



is expressed as [3 



Considering the quickly accumulating number of 7-ray 
sources observed by lACTs and Fermi, which will be im- 
proved essentially by the future Cerenkov telescope array 
(CTA) 27H29| , we propose to measure the EBL directly 
by extracting the collective absorption effects in a num- 
ber of 7-ray sources at different redshifts. In this work 
we demonstrate the capability of this method by adopt- 
ing the Markov Chain Monte Carlo (MCMC) fitting to 
known data to extract the parameters of the intrinsic 
spectra of a series of TeV blazars. as well as the EBL 
intensities simultaneously. Different from the previous 
studies in the literature, we make no assumption of the 
EBL model in our fitting. Instead the EBL is divided 
into many discrete energy bins and the energy density in 
each bin is fitted as a free parameter ^i. 

The observed VHE 7-ray spectrum after absorption by 
the EBL is commonly expressed as 



(1) 



where Fint{E) is the intrinsic spectrum of the source at 
redshift z. The strength of the attenuation by EBL is 
described by the optical depth t{E, z) as a function of 
energy E and the source redshift z. The optical depth r 



de'n'{e',z')a{E\e',fi), (2) 



where variables with prime are the quantities at redshift 



dl 



cdt 



dz' 



is the diffcren- 



"0 (1 + z')\/^m(1+z')-''+S^a 

tial path traveled by the VHE photon, fi ~ cos 6 with 9 
the angle between the momenta of VHE and EBL pho- 
tons, n'(e', z') = n{e'/{l + z'),z = 0)(1 + z')^ is the EBL 
number density at redshift 0', and a is the pair produc- 
tion cross section, e'^^ir is the threshold energy for 7-ray 
energy E' = E{1 + z') with an angle cos 6' = /i with the 
EBL photon. The cross section is peaked at a wavelength 
A//im ~ 1.24i;/TeV Therefore the observation of 

VHE 7-ray spectra can probe EBL at the wavelength 
from optical to far infrared, while it is not sensitive to UV 
band. The cosmological parameters used in this work are 
= 0.274, f^A = 1 - r^A/, Hq = 70.5 km s'^ Mpc'^ 

In the fitting the intrinsic spectra of blazars -Flnt are pa- 
rameterized by a power-law [F oc E~°') or log-parabolic 
{F cx E~°'~^^°s^) function, with two or three free pa- 
rameters for each source. The blazar spectrum is usually 
explained by the synchrotron-self-Compton (SSC) model, 
which shows a concave 7-ray spectrum. Such a spectrum 
can be formulated by the log-parabolic function. If the 
measured energy range is not too wide the simple power- 
law can give a quite good description. In the following 
we will compare the results with the two spectral forms. 

No prior assumption about EBL shape is adopted in 
this study. The EBL intensities are divided into 10 bins 
logarithmically between 0.1 and 100 fim. Within each bin 
the intensity vl^ is assumed to be a constant Then 
we can fit the 10 ^^s, according to Eq. ([T]), from a set of 
observed 7-ray spectra Fohs{E). 

We have adopted seven blazars in this study, which 
have relatively precise spectral measurements. The 
sources adopted are listed in Table ID Since the redshifts 
of these sources are less than 0.2, we neglect the evolution 
of the EBL 0. 

We employ the MCMC algorithm to do this global 
fit, which is very efficient for the minimization in high- 
dimensional parameter space [4o| . Physical constraints 
on the parameters are adopted. In models of diffusive 
shock acceleration of electrons in the blazar jets the de- 
duced 7-ray spectra are stro ngly constrained with power 
law index a larger than 1.5 [41[. The parameter (3 is re- 
stricted in [0, 1], and we shall test in the following that 
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TABLE I: Source sample information 
Name redshift Experiment Reference 



Mkn 421 


0.031 


VERITAS 


[33] 


Mkn 501 


0.034 


HEGRA 


[34] 


lES 1959+650 


0.047 


HEGRA 


[35] 


PKS 2005-489 


0.071 


HESS 


[36] 


PKS 2155-304 


0.116 


HESS 


[37] 


H 2356-309 


0.165 


HESS 


[38] 


lES 1101-232 


0.186 


HESS 


[16] 
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FIG. 1: The best-fitting results of the EBL intensities in 10 
wavelength bins. The intrinsic spectra of blazars are assumed 
to be log-parabolic and the low state of Mkn 421 spectrum 
by VERITAS is adopted. The solid lines show the upper 
limits derived by the 7-ray observations of blazars, HESS06 
[lil, Mazin07 [ig|, MAGIC08 [H, and Meyerl2 [3, and the 
dashed lines show the model predictions of Kneiske06 [jif . 
FranceschiniOS @|, FinkelO 01 and Gilmorel2 d. 



relaxing the range of does not change the results sig- 
nificantly. The EBL intensities are restricted in the 
lower and upper limits set by direct measurements 

Fig. [1] shows the fitting result assuming log-parabolic 
intrinsic spectra of balzars. There are seven states of Mkn 
421 observed by VERITAS [Hil, and the low state spec- 
trum (sample 2 as defined below) is adopted in this fit 
due to its wide energy coverage. For comparison, several 
recent model predictions of the EBL at 2; = Q-flEll are 
plotted by the dashed lines in the figure. Our fitting re- 
sult agrees very well with the model predictions for A > 1 
/im. The expected peaks at ~ 1/im and ~ lOOyum are 
also shown in this model-independent fit. It is very in- 
teresting that the indirect measurement of EBL provides 
strong support to the present understanding of EBL ori- 
gin. The solid lines in the same plot show several current 
upper bounds on the EBL intensities according to the 
7-ray observations of blazars 16l4l8l. |24 1 . 



In the following we discuss the robustness of the 
method. The largest uncertainty comes from the un- 
known nature of the source emission. Fortunately Mkn 
421 provides a perfect template to test the method. The 
observations by VERITAS of Mkn 421 covered seven 
states from very low to very high states from 2006 to 
2008 [3^. We can test the convergence of EBL by the 
fitting procedure using the different states data at the 
same source. Fig. [5] shows the results by combining the 
other 6 sources with the different states of Mkn 421 (re- 
ferred as sample 1-7 from very low to very high state) , for 
assumption of power-law (left) and log-parabolic (right) 
intrinsic spectra of the sources. For power-law intrin- 
sic spectra, the results show moderate diversity among 
different data set, while for the log-parabolic intrinsic 
spectra the results converge quite well. 

The reduced Xr = X^/d-O.f. are listed in Table HIl It 
is shown that for power-law source spectra the reduced 
varies for different data samples. In most of these 
cases the x^ values are too large. For the log-parabolic 
source spectra, the reduced Xr ~ 1 fo'^ almost all of the 
data samples. This means the SSC model predicted log- 
parabolic function can give a quite well description to 
the intrinsic 7-ray spectra. For Mkn 421, some of the 
seven observations span the energy range for more than 
1.5 orders of magnitude [33| . thus the simple power-law 
is not a good description. 

Then we relax our constraints on the parameter space. 
The EBL intensities are relaxed to be within [1,100] 
nW m^'^ sr^^. The spectral parameters are relaxed as 
a>2/3,0</3<2 respectively. The results are shown 
in Fig. [31 We notice that the current adopted data sam- 
ple can not constrain the EBL with A < 1 /J,m. The 
UV-optical band can be constrained by sources at high 
redshift and lower energy, such as the 3C 279 at z = 0.536 
17| . However, since only 5 points of 3C 279 was given 



by MAGIC it has no help to improve the fit because to 
include it we need to introduce 3 additional source pa- 
rameters. 

In summary we propose to measure the EBL from the 
VHE 7-ray data by a global fitting method. Both the 
intrinsic spectral parameters and the EBL intensities are 
fitted simultaneously using an MCMC algorithm, with- 
out any assumption of the spectral shape of EBL. With 
a log-parabolic VHE 7-ray spectra the fitting shows well 
convergence for the EBL intensities. The EBL intensities 
are close to the lower bound of EBL set by galaxy counts 
and are consistent with the recent EBL models. With 
the greatly improved number of VHE 7-ray sources by 
the future CTA the EBL can be determined with much 
higher precision. 
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FIG. 2: The best-fitting results of tlie EBL intensities in 10 wavelengtli bins, for tlie assumptions of power-law (left) and 
log-parabolic (right) intrinsic spectra of the sources. Different symbols are for different data sets (see the text for details). 

TABLE II: Best-fit x^/d-o.f. values 

1 2 3 4 5 6 7 

power-law 90.8/55 89.4/58 73.6/59 87.6/56 64.7/50 46.8/50 63.4/53 
log-parabolic 54.7/48 60.2/51 61.8/52 53.1/49 61.2/43 46.1/43 52.9/46 
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FIG. 3: Comparison of the results for relaxing fitting ranges of 
the EBL intensities or the spectral parameters. The reference 
case is the log-parabolic fit to data sample 2. 
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